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O ' We present a unique parameterization of the equation of state of strongly 

■ interacting matter in the temperature interval O.GTc ■ ■ ■ 3Tc at /x = within 

^ . a quasi-particle model based on quark and gluon degrees of freedom. The 

^ , extension to non-vanishing baryon-chemical potential is discussed. 

1 Introduction 

^ ■ The equation of state (EoS) of strongly interacting matter contains important in- 

CN . formation about bulk properties and the thermodynamics. In fact, one needs the 



X 



EoS as input for describing hydrodynamically the evolution of the early universe. 



in 

(N 

CN I dynamics and static characteristics of neutron stars as well as certain stages of rel- 

^ I ativistic heavy-ion collisions. Considerable progress has been achieved during the 

O ' last years in calculating the EoS from first principles resting on quantum chromo- 

^ . dynamics (QCD). By advanced sampling techniques the EoS of strongly interacting 

Q-i. matter is numerically accessible from lattice QCD calculations. Many details have 

been considered for particle-antiparticle symmetric matter. Here, at a certain value 
of the temperature, Tc, the susceptibilities corresponding to the Polyakov loop and 
the chiral condensate display pronounced peaks indicating deconfinement and chiral 
symmetry restoration (cf. [1 ). Above in the deconfinement region, the relevant 
degrees of freedom of strongly interacting matter are thought to be quarks and glu- 
ed ■ ons. Despite asymptotic freedom, in the range of < T < 2.5Tc, the EoS behaves 
highly non-trivial, preventing a direct perturbative treatment. There are various 
attempts to understand the EoS in this range and above, such as dimensional re- 
duction, resummed HTL scheme, $ functional approach, Polyakov loop model etc. 
(cf. [2j for a recent survey). A controlled chain of approximations from full QCD to 
some analytical expressions describing the lattice data without adjustable parame- 
ters would be of desire. Success has been achieved for T > 2.5Tc jSj- In contrast, 
the range T > Tc, in particular near to Tc, is covered by phenomenological models 
jU El El with parameters adjusted to lattice QCD data. 

Below Tc, in the confinement region, hadrons are thought to be the relevant 
degrees of freedom. Detailed information about the EoS in this range from lattice 
QCD became available very recently jH], offering a chance to check phenomenological 
models for the first time. In jU] lattice QCD data from [5| for T < Tc have been 
shown to agree with results of a hadron resonance gas model including a large 
number of states up to 2 GeV. Taking a different point of view, we present here a 
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description of the lattice QCD data [S] in the range 0.6Tc < T < 3Tc by our quasi- 
particle model of quarks and gluons |31 13] . Up to now it has been proven that our 
model can reproduce lattice data in the range Tc ■ ■ ■ ST^ for the pure gluon plasma 
[4j and for quark-gluon plasmas with various quark- flavor numbers 0. Here, our 
focus is to cover the lattice QCD data in the confinement region. 

Our paper is organized as follows. In section 2 we recapitulate features of our 
quasi-particle model. The comparison with lattice QCD data at vanishing chemical 
potential is performed in section 3. We discuss the obtained parameterization of 
the equation of state in section 4. Section 5 is devoted to non-vanishing chemical 
potential. The summary can be found in section 6. 



2 Quasi-particle Model 

Our quasi-particle model of light and strange quarks (g, s) and gluons (g) is based 
on the entropy density 



27r2T 



+ mf) 



mf + k"^ 



iUik) + f_{k)) - fiiuik) - f_{k)) 



(1) 



and the quark number density 

dq 
27r^ 

with degeneracies (i„ = 12, rf^ = 6 and rfg = 8 as for free partons, and distribution 



n = -JL dkk'[Uik)-f4k)] (2) 



functions f±{k) = {exp{[ymf + k'^ ^ lAI"^) + '^)~^ with S" = ( — 1) for quarks 
(gluons). The chemical potential is /i for light quarks and anti-quarks, respectively, 
while /i = for strange quarks and gluons. The parton effective masses 

T^1{T, /i) = mjo + Yii{k] T, /x) (3) 

have a rest mass contribution m?Q and, as essential part, the one-loop self-energies 
at hard momenta Ili{k;T, fi). The crucial point is to replace in Ili{k]T, fi) the 
running coupling by an effective coupling, G'^(T,fi) j4j. (As shown in p!U|, it is 
the introduced modification of G'^{T,fi) which allows to describe the lattice QCD 
data, while the use of the 1-loop or 2-loop perturbative coupling together with a 
more complete description of the plasmon term and Landau damping restricts the 
approach to T > 2.5Tc.) In doing so, non-perturbative effects are thought to be 
accommodated in this effective coupling. In massless $^ theory such a structure 
of the entropy density emerges by resumming the super-daisy diagrams in tadpole 
topology [TT], and [12] argues that such an ansatz is also valid for QCD. [3] point 
to more complex structures, but we find (^Q ^ flexible enough to accommodate the 
lattice data (see below). 

The pressure reads accordingly 

p= ^ - i?(T,/x), P. = ^ / dk [Uik) + f_ik)] , (4) 
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where B{T,fi) ensures thermodynamic self-consistency, s = dp/dT, n = dp/dfi 
together with the stationarity condition Sp/6m] = [TH]. Note that eqs. ([TJ |21 
themselves are highly non-perturbative expressions: Expanding (jU in powers 
of the coupling strength one recovers the first perturbative terms. (Higher order 
terms would probably need higher orders in 11). For more details on our model see 

H 001101112]. 



3 Comparison with lattice data at /i = 

represents a mapping of lattice data for s(T) on G^(T); to determine p one has 
to fix an integration constant, say B{T^. As in the lattice calculations [8^ we have 
used mj = xiT with Xq = 0.4 for light quark flavors, Xg = 1 for strange quarks and 
Xg = 0. We have found as a convenient parameterization 

G^(T) = P^-'-^' (5) 
[ a - bT/T„ T<% 

where Gg^ioop is the 2-loop coupling 

2 ^ 167r^ _ 



2Alog(loga 



(6) 



with (3q = {IIN^ - 2Nf)/3, (3i = {SAN^ - UNfN^ + 3Nf/Nc)/6, and the argument 
^ = A(T — Ts)/Tc. The parameters are Tg = O.SOTc, A = 7.6, a = 359.0 and 
b = 334.5 for Nf = 2 + 1, Nc = 3. The comparison of our model with the lattice 
data |HJ inj is exhibited in Fig. 1 by solid curves for entropy density and pressure, 
respectively. One observes an impressively good description of the lattice data, in 
particular also below Tc. 

Remarkable is the change of the shape of G'^{T) at Tc, see Fig. 2. Above Tc, 
G'^{T) resembles the usual 2-loop running coupling strength with a regulator T^, 
which ensures agreement with the data down to Tc. Going down in T, at Tc the 
growing of G'^{T) is changed into a moderate linear rise. Note that no order pa- 
rameter is needed which changes at Tc. This is possibly related to the fact that s 
is a measure for the density of states. Also, the degeneracies di are constant. 

Here we would like to contrast our quasi-particle model based on quark-gluon 
degrees of freedom with the resonance gas model below Tc jS]. The authors of jH] 
state that many heavy resonances are needed to describe the data (at Tc the lightest 
hadron contribution to the energy density is at the level of 15% [H]). Our result 
is in some correspondence, in the sense that fairly heavy quasi-particle excitations 
are needed which emerge from the large values of G^. One could also argue that a 
large number of excited hadron states at T < Tc can be effectively described in an 
equivalent manner by a small number of quasi-particle excitations, as we do here. 

It is well conceivable that, for determining the pressure (j3)), two distinct inte- 
gration constants would be needed, say B(T^) and B{T^). However, one constant 
B{Tc) = OAAT^ is sufficient to describe the pressure, see right panel in Fig. [TJ 
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As mentioned above, the lattice calculations [H| use quite heavy quarks. A possi- 
ble estimate of the chiral extrapolation is to put rriifi = and repeat the evaluation 
of eqs. (P-IH) with keeping the parameters in (jSJ fixed; B(Tc) is adapted to en- 
sure a positive pressure. The results are shown as dashed curves in Fig. ^ Such a 
procedure is certainly too simple. Only once lattice results for various rriifi are at 
our disposal one can estimate the mj q dependence of the parameters in (jSJ and for 
B(Tc) to achieve a more profound chiral extrapolation. 
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Figure 1: Entropy density s scaled by (left panel) and pressure scaled by (right 
panel) as a function of the temperature in units of the pseudo-critical temperature 
Tc- The full (dashed) curves represent calculations with lattice masses as described 
in text (rrtifl = as an estimate of the chiral extrapolation). Lattice data from 
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Figure 2: The effective coupling G'^{T) (in units of G'^{Tc)) as a function of the 
temperature. 



4 Discussion 



Up to this point we consider eqs. (P - 0} as a convenient unique parameterization 
of the EoS, which one may make use of for purposes mentioned in the introduc- 
tion. Further detailed lattice data are needed to test our model at even smaller 
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temperatures and to see where an explicit change to hadronic degrees of freedom is 
required. 

In an restricted interval around one may speculate whether a quark-hadron 
duality is at work, i.e. hadron observables are expressed in a quark-gluon basis and 
vice versa. For example, jU] shows that in a narrow interval above the hadron 
resonance gas model still agrees with the lattice data. Above T^, signals of hadronic 
states have been found as well ^3] , an aspect advocated as important in ^3] . Even 
more, [IB] has shown that a modified resonance gas model describes the lattice data 
above Tc also successfully. 

There are various examples of the quark-hadron duality in the literature. We 
mention here QCD sum rules [El which realize the duality, as discussed in [TH] . 
Another observation is that in heavy-ion collisions the low-mass di-electron spectra 
can be described by a quark-gluon plasma emission rate, even if the space-time 
averaged temperature of the emitting system is below ^j. For further aspects 
of the duality we refer, e.g., to [2Uj . 

5 Non-vanishing chemical potential 

With the advent of lattice QCD data at non- vanishing chemical potential [2111221123] 
the equation of state becomes accessible in a large region. This is particularly impor- 
tant as the envisaged CBM project at the future accelerator complex in Darmstadt 
aims at exploring systematically the region of maximum baryon densities at 
reach in heavy-ion collisions. 

As described in jHl [12] our quasi-particle model also covers the equation of state 
for non-vanishing chemical potential. The point here is that thermodynamical self- 
consistency allows to map the lattice QCD data from /i = to > without further 
assumptions. Indeed, as shown in |2S] our model agrees quite perfectly with lattice 
QCD calculations at yU > for 2 + 1 flavors. Here we focus on the recent lattice QCD 
data for 2 flavors [2S1 with improved p4 action. Given G'^(T) (e.g., from the previous 
section, or from an analysis of other lattice QCD data such as |2S|) as boundary 
values, the partial differential equation ardG'^iT, jj) / dT + a^dG^iT, fi)/dfi + aT,^j. = 
with aT(/^ = 0) = and a^(T = 0) = (cf. [HI [12] for details) delivers as 
solution G^{T, /x), needed in Ili{k; T, /x). Fig. 3 exhibits the characteristic curves for 
solving this "flow equation" . Asymptotically, G^ is constant along the characteristic 
curves. As discussed in [HI E] we take the characteristic curve emerging from 
a.t n = (solid curve in Fig. 3) as indicator of the phase border line. It agrees 
up to intermediate values of n fairly well with estimates based on lattice QCD 
calculations [22.. The characteristic curves above the solid curve show a regular 
pattern of sections of ellipses. In contrast, the characteristic curves emerging from 
the T axis below are flatter and cross the former ones at larger values of fi. This 
may indicate that the model is not longer applicable in that region. To understand 
the pattern of the characteristic curves we note that, in the region of interest, the 
characteristic curves have a universal shape, i.e., for small values of fi the shape is 
independent of the actual start value G'^{T,^ = 0). For sufficiently small values of 
G'^{T,fi = 0), the characteristic curves form a nested set of non- intersecting ellipse 
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like curves. Increasing G^{T,fi = 0) causes a levelling-off at larger values of fi. The 
levelling-off sets in at smaller values of fi for larger values of G'^{T,fi = 0). This 
causes the pattern seen in Fig. 3. 

Having G^{T,fj,) at our disposal, the EoS follows from eqs. (1 - 4) in a straight- 
forward way. Fig. 4 shows as example a comparison of our model with the lattice 
QCD data from for Nf = 2? Again, a fairly well reproduction of the data 
is achieved for T > Tc, while at T < T^, in particular for large values of /i, the 
agreement is less perfect. 




Figure 3: The characteristic curves solving the flow equation for the effective cou- 
pling . Solid (dashed) lines: characteristic curves emerging from T > (T < 
Tc). The heavy solid curve is considered as indicator of the phase border line. The 
characteristic curves are not drawn in the region where they intersect. The marked 
region indicates where the lattice QCD data are given. Note that they are ac- 
tually provided by the truncated expansion Ap/T'^ = p{T, fi = 0)/T^ —p(T, fJ')/T^ = 
C2{T){^/TY + Cis^{T){jj,/TY . Higher order terms, like Cq(T) (fi/T)^ etc. are needed 
to control the validity of the extrapolation to larger values of fi. 



6 Conclusions 

To summarize, we present a unique parameterization of the equation of state of 
strongly interacting matter in the temperature interval O.GTc ■ ■ ■ 3Tc which is based 
on the picture of weakly interacting quasi-particles of quarks and gluons. At Tc 
for = 0, the temperature dependence of the effective coupling suffers a change; 
no rapidly changing order parameter is needed to describe the lattice QCD data. 
Further precision lattice data below are needed to check the model in more 
detail and to put the speculation on firm ground, such that the applicability of the 
quasi-particle picture provides an additional example of quark-hadron duality. 

^Here the parameters for the effective couphng in eq. (5) are given by A = 21.3, Tg — 0.93Tc, 
a = 427.4, and b = 404.2. 
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Figure 4: The scaled pressure difference Ap(T, jj) = p(T, jj) —p(T, /i = 0) (left panel) 
and the scaled quark number density (right panel) as a function of the temperature 
for various chemical potentials. Data from \2!^ . 



It has been shown that our model is also able to describe the new lattice QCD 
data for non- vanishing chemical potential, as alternative approaches do pm i27j. As 
particular new point we note that we cover the region T < at non-vanishing but 
small chemical potential with our quasi-particle model. 

While present lattice QCD data deliver either p{T,fi = 0) or Ap(T, /i), our 
model covers both quantities on the same footing. 

Inspiring discussions with J. P. Blaizot, D. Blaschke, A. Peshier, and K. Redlich 
are gratefully acknowledged. The work is supported by BMBF 06DR121, GSI and 
FP6-I3 network. 
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